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ABSTRACT 

When  a  circular  cylinder  is  immersed  in  a  fluid  of  infinite  extent 
and  oscillated  in  such  a  manner  that  drag  effects  are  negligible,  it  is 
well  known  that  the  closed  form  solution  gives  an  added  mass  coefficient 
of  1.0.  However,  when  boundaries  are  placed  in  proximity  to  the  cylinder, 
and  the  cylinder  oscillated  parallel  to  the  boundary,  the  added  mass 
varies  from  the  classical  result.  The  forces  acting  on  the  cylinder  in 
the  latter  case  are  dependent  upon  the  proximity  of  the  boundaries  as 
well  as  the  frequency  of  oscillation. 

It  was  the  purpose  of  this  study  to  experimentally  determine  the 
total  force  on  such  a  cylinder  in  proximity  to  an  impermeable  bottom 
surface  and  under  a  free  surface  and  to  relate  this  force  to  the  added 
mass  coefficient.  Experimental  data  was  reduced  using  a  computer  program 
and  the  results  presented  graphically  as  a  function  of  dijiaensionless 
parameters . 
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NQMENCLATUKE 
Symbol  Description 

a  Radius  of  Circular  Cylinder 

d  Relative  Depth  of  Submergence 

3"  Depth  of  Submergence  of  Cylinder 

F  Total  Maximum  Force  Per  Unit  Length 

g  Acceleration  of  Gravity 

h  Relative  Water  Depth 

h  Water  Depth 

PA  Phase  Angle 

WH  Waveheight 

x  Relative  Displacement 

X  Half  Amplitude  of  Cylinder  Oscillation 

v  Dimensionless  Frequency  Parameter 

p  Water  Density 

a  Circular  Frequency  of  Cylinder  Oscillation 

y  Water  Viscosity 


DIMENSIONLESS  NOTATION 

F/pgTra2cr2X  Force  Coefficient 

a2a/g  =  v  Frequency  Parameter 

d/a  =  d  Relative  Depth  of  Submergence 

h/a  =  h  Relative  Water  Depth 

X/a  =  x  Relative  Displacement 

2paXcr/y  Reynold's  Number 
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I.   INTRODUCTION 

During  the  past  two  decades,  many  new  submarine  structures,  such  as 
pipelines,  oil  storage  tanks  and  large  aquaducts  have  been  either  sunk 
or  moored  on  the  continental  shelf.  Consequently,  the  calculation  of  the 
forces  exerted  on  such  objects  is  of  practical  importance. 

In  an  early  paper  [1] ,  Stokes  showed  that  the  expression  for  the 
force  on  a  pendulum-like  system  oscillating  in  an  infinite  fluid  was  both 
acceleration  and  velocity  dependent. 

In  more  recent  years,  Morison,  et  al,  [2] ,  from  an  investigation  of 
wave  forces  exerted  on  piles,  showed  that  when  the  characteristic  dimen- 
sion of  a  body ,  i.e.,  pile  diameter ,  was  small  compared  to  wave  length , 
the  force  consisted  of  both  a  drag  and  an  inertial  term.  However,  the 
coefficient  of  each  of  these  terms  depends  upon  the  geometry  of  the 
object  in  question,  the  amplitude  of  the  fluid  motion,  the  proximity  of 
a  free  surface,  and  the  proximity  of  an  impermeable  bottom  surface. 

Several  investigations  have  been  conducted  since  Morison's  work,  which 
have  dealt  with  sinusoidal  flow  about  various  objects.  Keulegan  and 
Carpenter  [3]  performed  an  experimental  study  with  cylinders  and  plates 
in  a  rectangular  tank  in  which  a  standing  wave  was  generated  in  water. 
Forces  on  the  submerged  bodies  were  measured  with  a  dynamometer.  Inertial 
and  drag  coefficients  were  correlated  with  a  dimensionless  "period 
parameter"  for  both  unseparated  and  separated  flow,  the  period  parameter 
being  proportional  to  the  relative  displacement  of  the  surrounding  fluid. 
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Heinzer  [4]  carried  out  a  visual  study  of  the  flow  field  character- 
istics of  a  cylinder  sinusoidally  oscillated  in  a  still  fluid.  Flow 
patterns  included  laminar  boundary  layer,  symmetric  separation,  and 
Karman  vortex  street.  Heinzer 's  major  conclusion  was  that  wake  charac- 
teristics were  primarily  dependent  upon  the  relative  amplitude  of  oscil- 
lates. A  consideration  of  this  type  was  used  in  the  study  which  follows. 
Snail  amplitude  oscillations  were  utilized  at  various  frequencies  to 
insure  unseparated  flow. 

Sarpkaya  and  Garrison  [5]  studied  constantly  accelerated, 
unidirectional  flow  past  a  circular  cylinder.  Forces  were  measured 
experimentally  and  the  nature  of  the  flow  determined  visually.  Force 
coefficients  were  calculated  as  functions  of  the  relative  displacement 
of  the  cylinder.  It  was  shown  that  the  total  resistance  as  well  as  the 
drag  and  inertial  forces  are  representable  as  functions  of  the  relative 
displacement  of  the  fluid. 

Ward  and  Dalton  [6]  analyzed  the  forces  exerted  on  a  cylinder  in  a 
purely  sinusoidal  flow  field  in  which  symmetrically  located  vortex  pairs 
were  developed.  Forces  acting  on  the  cylinder  were  calculated  using 
Bryson's  [7]  potential  flow  model  including  a  pair  of  vortices  but  no 
attempt  was  made  to  verify  these  forces  experimentally. 

Later  Dalton  and  Hamann  [8]  conducted  an  experimental  investigation 
of  a  sinusoidally  osciallating  cylinder  in  a  quiescent  fluid  in  an  attempt 
to  represent  a  sinusoidal  wave  moving  past  a  single,  fixed  platform  leg. 
The  conclusions  drawn  from  their  investigation  were  mostly  qualitative. 
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The  purpose  of  this  study  was  to  provide  a  correlation  between 
the  added  mass  coefficient  and  appropriate  diirensionless  parameters.  A 
right  circular  cylinder  was  oscillated  in  a  quiescent  fluid  at  ampli- 
tudes which  insured  unseparated  flow.  The  results  are  presented  in 
dimensionless  form  as  a  function  of  a  dimensionless  frequency  parameter, 
relative  water  depth,  relative  depth  of  submergence,  and  relative 
amplitude  of  oscillation. 
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II.   THEORETICAL  CONSIDERATIONS 

A  schematic  representation  of  the  problem  is  shown  in  Figure  (1) .  A 
right  circular  cylinder  of  radius  a  was  submerged  to  a  depth  d  in  water 
of  depth  h.  The  cylinder  was  oscillated  sinusoidally  at  a  circular  fre- 
quency a   with  half  amplitude  of  oscillation  x   Of  primary  interest  was 
the  total  maximum  force  per  unit  length  F  exerted  by  the  water  on  the 
cylinder,  the  phase  of  this  force  and  the  amplitude  of  the  generated  waves. 

The  complete  theoretical  solution  to  the  problem  is  complicated,  even 
when  the  amplitude  of  oscillation  is  kept  small  enough  to  justify  an 
assumption  that  the  flow  remained  unseparated.  Therefore,  in  order  to 
make  seme  headway  a  dimensional  analysis  was  applied  to  the  problem.  It 
is  known  "a  priori"  that  the  maximum  force  per  unit  length  of  cylinder  has 
the  following  functional  dependence. 

F  =  £^(h,d,a,a  ,X',p  ,qf\i)  (1) 

where  h  denotes  the  water  depth,  d  the  depth  of  submergence,  a  the  cylinder 
radius,  a   the  circular  frequency  of  the  imposed  motion,  X  the  half  ampli- 
tude of  the  imposed  motion,  p  the  water  density,  g  the  acceleration  of 
gravity,  and  y  the  water  viscosity. 

Dimensional  analysis  of  the  physical  quantities  of  Equation  (1) 
resulted  in  the  following  dimensionless  groups 

r./   -?  2v   ^  ,a2a  d  h  X   2paX?\  ,0* 

F/pgiraAjOC  =  f0(— —  /  =,  -,  =,    — - — )  (2) 

^  9       a  a  a    v 
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x  =  Xsin(Tt 
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FIGURE  1.   DEFINITION  SKETCH 
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«     2paXj  .   , 
where  ~ is  the  well  known  Reynold '  s  number  representing  the  ratio  of 

inertial  to  viscuous  forces.  Since  the  cylinder  was  oscillated  at  small 
amplitudes  it  was  assumed  that  flow  around  the  cylinder  was  unseparated 
and  viscous  forces  small.  Therefore/  the  Reynold's  number  was  disregarded 
and  Equation  (2)  written  as 

F/pgriPa*  -  f,.^  £!L:§  (3) 

J  y   a  a  a 

A  second  reason  for  ignoring  a  dependence  upon  Reynold's  number  was 
provided  by  Heinzer  [4]  .  His  visual  studies  showed  that  wake  characteris- 
tics were  primarily  dependent  upon  the  relative  amplitude  of  oscillation 
and  that  Reynold's  number  effects  were  not  evident. 

As  long  as  the  relative  displacement  is  small,  the  flow  remains 
essentially  unseparated  and,  accordingly,  the  force  is  also  independent 
of  the  relative  amplitude.  The  independence  of  the  force  from  this  para- 
meter for  small  values  was  demonstrated  experimentally  by  Keulegan  and 
Carpenter  [3]  .  In  the  present  investigation,  no  attempt  was  made  to  study 

the  effects  of  this  parameter;  consequently,  it  was  held  fixed  at  —  =  0.1 

a 

throughout  the  experiment. 

The  objective  of  this  experimental  work  was  to  show  a  correlation 
between  the  dimensionless  force  coefficient  on  the  left  hand  side  of 
Equation  (3)  and  the  dimensionless  parameters  on  the  right  hand  side. 
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III.   EXPERIMENTAL  INVESTIGATION 

A.  APPARATUS 

1.  Wave  Channel 

The  wave  channel  was  designed  and  built  by  Shiller  [9]  for  an 
earlier  thesis.  The  channel  was  designed  with  a  rectangular  cross  section 
two  feet  deep  and  fifteen  inches  wide.  It  was  fabricated  in  sections, 
each  section  consisting  of  two  vertical  walls  and  one  horizontal  floor 
member.  The  vertical  walls  were  stiffened  by  use  of  three  four  foot  long 
2  x  4's  attached  to  3/4  inch  x  8  foot  x  3  foot  sheets  of  plywood  at  the 
ends  and  middle  of  the  longest  dimension.  The  channel  floor  member  was 
composed  of  two  15  inch  x  8  foot  long  plywood  sheets  separated  by  a 
1  inch  x  4  inch  strip  around  the  edge.  The  bottom  and  side  members  for- 
formed  the  basic  channel  cross  section.  The  bottom  and  side  members  were 
mated  together  with  3/8  inch  threaded  rods.  Six  of  these  8  foot  sections 
were  sealed  and  bolted  together  to  form  a  48  foot  long  channel.  Square- 
ness of  the  channel  cross  sections  was  maintained  by  use  of  spacers  placed 
between  the  stiffened  members  near  the  floor  and  across  the  top. 

One  end  of  the  channel  was  fitted  with  a  wave  generating  mechanism 
which  was  not  used  in  this  investigation. 

At  a  location  30  feet  from  the  wave  generator  end  of  the  channel, 
a  three  foot  wide  plexiglass  observation  window  extending  the  full  height 
of  the  channel  was  installed  in  both  walls  as  shown  in  Figure  (2) .  The 
test  cylinder  was  sinusoidally  oscillated  at  a  position  near  one  end  of 
the  plexiglass  window. 
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The  generated  waves  were  dissipated  on  a  beach  structure  located 
at  each  end  of  the  channel.  The  beach  as  shown  in  Figure  (3)  consisted 
of  three  layers  of  perforated  stainless  steal  sheeting.  Layer  separation 
and  structural  strength  were  provided  by  2  x  4  runners  with  their  largest 
cross  sectional  dimension  perpendicular  to  the  sheeting.  The  slope  of 
the  beach  was  6:1. 

2.  Carriage  and  Driving  Mechanism 

The  driving  mechanism  as  shown  in  Figure  (4)  consisted  of  a  vari- 
able speed  transmission  driven  by  a  two-horsepower  inductance  motor.  The 
motor  and  transmission  were  mounted  on  a  3/4  inch  steel  plate  which  was 
positioned  and  bolted  to  the  top  of  the  channel.  The  plate  was  further 
braced  using  1/2  inch  angle  iron  supports  which  were  bolted  across  the 
top  of  the  channel  and  diagonally  to  the  channel  side  wall.  The  variable 
speed  transmission  provided  an  output  shaft  speed  range  extending  from  22 
to  180  rpm  which  in  turn  provided  dijrensionless  frequency  parameters  of 
from  about  0.04  with  a  six  inch  diameter  cylinder  to  about  1.7  with  a 
four  inch  diameter  cylinder. 

The  output  shaft  of  the  transmission  was  fitted  with  a  face  plate. 
The  face  plate  as  shown  in  Figure  (5)  was  connected  to  the  drive  rod 
through  a  bolt  and  slot  assembly  which  provider1,  a  drive  rod  stroke  which 
could  be  varied  between  zero  and  six  inches. 

The  carriage  as  shown  in  Figure  (6)  provided  a  movable  mounting 
platform  capable  of  being  positioned  over  the  top  of  the  wave  channel. 
The  carriage  was  driven  sinusoidally  by  the  driving  mechanism  at 
amplitudes  of  about  10%  of  the  cylinder  diameter. 
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FIGURE  3.   BEACH  STRUCTURE 
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FIGURE  5.   FACE  PLATE  WITH  BOLT  AND  SLOT  ASSEMBLY 
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The  carriage  consisted  of  a  U-shaped,  welded  aluminum  platform 
fitted  with  four  3/4  inch  linear  bearings.  The  sides  of  the  platform 
were  constructed  of  two  inch  OD  aluminum  pipe  fitted  at  each  end  with 
aluminum  bearing  blocks  designed  to  accommodate  the  linear  bearings. 
The  structural  transverse  member  of  the  carriage  was  a  2  x  3  inch  alumi- 
num box  which  was  welded  to  the  pipes  at  one  end.  The  cantilever  beam 
used  to  measure  forces  on  the  cylinder  was  mounted  on  this  transverse 
member.  At  the  opposite  end  of  the  aluminum  pipe,  the  bearing  blocks 
were  connected  by  a  1/2  inch  diameter  transverse  shaft  which  accommodated 
the  cylinder  holding  and  balancing  jig. 

The  U-shaped  platform  was  mounted  through  the  linear  bearings 
onto  two  case-hardened  3/4  inch  shafts.  The  shafts  were  supported  and 
held  parallel  by  two  six  inch  channel  end  pieces  which  transversed  the 
wave  channel.  The  six  inch  channel  end  pieces  were  clamped  to  the  top  of 
the  wave  channel  at  a  position  over  the  observation  window.  The  driving 
mechanism  oscillated  the  carriage  sinusoidally.  The  driving  rod  was 
connected  to  the  2x3  inch  aluminum  transverse  member  by  a  stiff  link 
which  passed  through  an  accommodation  hole  drilled  through  the  intervening 
six  inch  channel  end  piece. 

3.  Cylinder  Holding  and  Balancing  Jig 

The  cylinder  holding  and  balancing  jig  consisted  of  two  long  arms 
pivoted  near  their  centers  as  shown  in  Figure  (7) .  The  jig  was  designed 
to  hold  the  test  cylinder  at  the  desired  depth  of  submergence.  It  was 
designed  to  have  as  low  a  moment  of  inertia  with  respect  to  its  pivot 
point  as  possible  in  order  to  maintain  a  high  system  natural  frequency 
while  having  adequate  sensitivity  of  force  measurement.  A  counter  balance 
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was  placed  on  the  arms  opposite  the  cylinder  and  adjusted  so  that  the 
center  of  gravity  of  the  system  was  at  the  pivot  point.  As  a  result  of 
this  balance  adjustment  the  only  force  (or  moment)  acting  on  the  arms 
was  that  due  to  the  water;  all  other  inertial  loads  were  balanced  out. 

A  twelve  inch  long  arm  was  welded  at  the  longitudinal  center  of 
the  center  drum  and  functioned  to  convert  a  torque  about  the  rotational 
center  of  the  drum  to  a  linear  force.  This  force  acted,  through  a  ball 
bearing  linkage,  on  a  small  cantilever  beam.  Two  attachment  points  on 
the  beam  were  provided  so  that  a  large  range  of  forces  could  be  measured 
without  overstressing  the  beam.  The  stiffness  of  the  beam  was  adjusted 
for  adequate  sensitivity  while  providing  the  maximum  possible  system 
natural  frequency.  The  design  natural  frequency  was  approximately  5.0 
tines  the  maximum  forcing  frequency. 

4.  Beam  and  Link 

The  beam  and  link  assembly  as  shown  in  Figure  (8)  functioned  to 
rigidly  connect  the  cylinder  holding  and  balancing  jig  to  the  carriage. 
Strain  gauges  were  applied  to  the  beam  to  provide  a  means  of  measuring 
forces  imposed  upon  the  cylinder  and  transmitted  mechanically  through 
the  assembly.  The  strain  gauge  output  was  calibrated  using  dead  weight 
loads  applied  horizontally  to  the  cylinder  through  a  pully  arrangement. 

5.  Test  Cylinders 

Two  cylinders  as  shown  in  Figure  (9)  were  constructed  to  provide 
maximum  flexibility  in  testing  over  a  wide  range  of  dimensionless  fre- 
quency parameters.  A  four  inch  CO  cylinder  was  constructed  of  thin- 
walled  aluminum  tubing  and  a  six  inch  OD  cylinder  was  constructed  of 
plexiglass.  The  cylinders  were  slotted  to  accommodate  one  end  of  the 
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FIGURE  9.      TEST  CYLINDERS 
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holding  and  balancing  arms.  Bulkheads  were  placed  adjacent  to  the  slots 
to  provide  an  attachment  point  for  the  ends  of  the  arms.  Silicone  sealant 
was  used  to  seal  the  cylinder  against  leakage  after  the  arms  were  attached. 
This  allowed  the  cylinder  to  be  unflooded  resulting  in  a  ndnimum  fixed 
mass. 

End  pieces  as  shown  in  Figure  (10)  were  constructed  of  plexiglass 
for  each  cylinder.  The  end  pieces  were  fitted  into  the  ends  of  the 
cylinders  with  1/8  inch  "0"  rings  to  provide  the  necessary  seal.  End 
pieces  were  "dished"  on  the  inner  surface  to  reduce  weight.  The  outer 
surfaces  of  the  end  pieces  had  a  modified  labyrinth  seal  cut  into  the  face 
to  prevent  water  from  flowing  between  the  wave  channel  wall  and  the 
cylinder  end  when  the  cylinder  was  oscillated. 
6.  Wave  Height  Gauge 

The  wave  height  gauge,  as  shown  in  Figure  (11) ,  consisted  of  a 
small  diameter  aluminum  rod  which  acted  as  the  basic  frame  to  which  a 
plexiglass  mounting  table  was  attached.  The  elements  of  a  Wheatstone 
bridge  were  mounted  on  this  table,  and  two  30  gauge  copper  wires  spaced 
1/2  inch  apart  were  extended  from  the  table  and  attached  to  an  insulator 
mounted  at  the  end  of  the  aluminum  rod.  The  other  end  of  the  rod  was 
fitted  to  a  block  and  stand  which  allowed  the  gauge  to  travel  ±6  inches 
in  the  vertical  direction.  The  block  and  stand  provided  a  convenient, 
accurate  means  of  calibration  through  adjustment  of  the  submergence  of 
the  wires.  The  complete  unit  was  positioned  over  the  channel  and 
centered  at  about  five  feet  from  the  carriage. 

The  wave  height  gauge  was  of  the  parallel  wire  resistance  type. 
The  gauge  operated  on  the  principle  that  the  conductance  between  two 
parallel  submerged  wires  varies  proportionately  to  the  length  of  wire 


28 


n 
w 


H 


29 


FIGURE  11.   WAVE  HEIGHT  GAUGE 
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submerged .  The  wave  height  gauge  was  connected  in  parallel  to  one  leg 
of  a  Wheatstone  bridge  as  shown  in  Figure  (12) .  To  conpensate  for  the 
inability  of  the  carrier  preamplifier  to  offset  the  resistance  of  the 
immersed  parallel  wires,  an  external  fixed  resistance  was  connected 
across  an  opposite  leg  of  the  bridge.  The  Wheatstone  bridge  was  connected 
with  shielded  cable  to  a  carrier  preamplifier  and  recorder. 

7.  R.  P.  M.  Switch 

The  rpm,  switch  as  shown  in  Figure  (13)  was  designed  to  measure 
the  frequency  of  the  motion  imposed  on  the  cylinder.  At  the  end  of  each 
stroke,  the  switch  was  closed  and  the  event  recorded  on  the  chart.  This 
mark  was  also  used  to  determine  the  phase  shift  between  the  force  and 
acceleration. 

8.  Carrier  Preamplifier  and  Recorder 

A  two-channel  Hewlett  Packard  8805B  carrier  preamplifier  and 
recorder  was  used.  One  channel  recorded  force  on  the  cylinder  and  the 
second  recorded  wave  height.  The  event  marker  was  utilized  through  an 
external  jack  to  record  the  frequency  of  the  imposed  motion  and  to 
establish  the  motion/force  phase  relationship. 

B.   TEST  PROCEDURE 

1.  System  Selection  and  Optimization 

The  first  step  in  the  procedure  was  to  select  one  of  two  sets  of 
holding  and  balancing  arms  to  provide  the  desired  d/a  ratio  and  to  mount 
them  on  the  center  drum.  Wooden  shims  were  available  in  1/2  inch  thick- 
ness which  allowed  placement  of  either  cylinder  at  two  or  three  radii 
from  the  wave  channel  bottom.  Once  a  set  of  holding  and  balancing  arms 
was  chosen,  the  dimensionless  parameter  d/a  varied  with  water  depth. 
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FIGURE  13.      R.P.M.    SWITCH 
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The  next  step  was  to  decide  whether  it  was  desired  to  examine  the 
lower  or  higher  range  of  dimensionless  frequency  parameters.  The  larger 
six  inch  OD  plexiglass  cylinder  was  used  for  the  lower  range  of  frequency 
parameters  because  its  larger  displaced  volume  produced  forces  which  were 
more  readily  measurable.  However,  at  higher  frequencies  the  larger 
cylinder  resulted  in  natural  to  forcing  frequency  ratios  too  close  to  1.0 
to  give  valid  results.  The  larger  six  inch  OD  cylinder  was,  therefore, 
used  at  the  lower  motor  speed  range  of  from  22  to  about  100  rpm.  The 
smaller  four  inch  OD  aluminum  cylinder  was  used  at  upper  speed  range  of 
from  about  100  to  180  rpm. 

Once  the  combination  of  holding  and  balancing  arms  and  cylinder 
was  chosen,  the  cylinder  was  attached  to  the  arms.  The  cylinder  was 
sealed  and  the  cylinder  holding  and  balancing  jig  was  balanced  using  the 
adjustable  counter  balancing  weight.  After  balancing,  the  jig  was  placed 
in  position  on  the  carriage  and  the  carriage  properly  positioned  and 
clamped  to  the  top  of  the  wave  channel. 

The  desired  stroke  was  adjusted  at  the  variable  speed  trans- 
mission face  plate.  The  stroke  was  maintained  at  a  value  of  0.1  a  to 
insure  that  separated  flow  did  not  occur. 

The  channel  was  filled  and  data  runs  conducted  at  water  levels 
which  were  2,  3,  and  4  cylinder  radii  above  the  cylinder  center. 

After  the  wave  channel  was  filled  to  the  desired  level,  the 

Wheatstone  bridges  were  balanced  and  a  test  run  was  made  through  the 

frequency  range  to  determine  the  optimum  attenuation  settings,  chart 

speed,  and  beam  position  to  obtain  the  most  readable  results.  At  this 

time  the  entire  assembly  was  "plucked"  to  determine  the  natural  to  forced 

frequency  ratio.  If  the  value  of  this  ratio  was  less  than  five,  the  beam 

was  adjusted  to  the  position  which  provided  the  stiff  est  beam. 
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2.  System  Calibration 

a.  Force  Calibration 

After  balancing  the  bridge  and  selecting  the  appropriate 
attenuator  setting,  a  piece  of  light  weight  plastic  line  was  tied  loosely 
around  the  center  of  the  test  cylinder  and  bent  around  a  low  friction 
pulley  located  at  the  top  of  the  channel,  100  inches  from  the  test  cylinder. 
Weights  were  hung  from  the  line  and  the  resulting  deflection  on  the 
recording  paper  recorded.  The  angle  which  the  string  made  with  the  hori- 
zontal was  accurately  measured  and  the  resultant  horizontal  component  of 
force  determined.  The  strain  gauge  bridge  was  calibrated  in  this  manner 
in  both  directions  separately.  The  results  were  plotted  accurately  to 
insure  linearity  and  to  be  used  to  interpret  data.  The  calibration  of 
the  strain  gauge  bridge  was  checked  frequently  between  data  runs.  See 
Appendix  A  for  a  typical  force  calibration  plot. 

b.  Wave  Height  Calibration 

After  balancing  the  bridge  and  selecting  the  proper  attenu- 
ator setting  the  wave  height  gauge  was  calibrated  by  utilizing  the  block 
stand  to  immerse  and  withdraw  the  wire  a  known  distance.  The  following 
procedure  was  used.  The  gauge  was  raised  1/2  inch  simulating  a  decrease 
in  water  level,  then  lowered  1/2  inch  below  the  null,  simulating  an 
increase  in  water  level.  Then,  the  gauge  was  raised  one  inch  and  lowered 
one  inch,  etc.  This  procedure  closely  simulated  the  actual  passage  of  a 
wave  profile  with  crests  and  troughs.  Resulting  deflections  on  the 
recording  paper  were  accordingly  recorded  and  the  results  accurately 
plotted  to  insure  linearity  and  to  be  later  used  to  interpret  data.  The 
calibration  of  the  wave  height  meter  was  checked  frequently  between  data 
runs.  See  Appendix  A  for  a  typical  wave  height  calibration  plot. 
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c.  Event  Marker  Calibration 

The  ability  of  the  event  marker  to  record  the  frequency  of 
the  imposed  motion  and  the  phase  angle  was  dependent  upon  the  accuracy 
of  the  chart  speed.  The  accuracy  of  the  chart  speed  was  checked  against 
the  internal  timer  of  the  recorder.  No  discrepancy  was  found. 
3.  Data  Run 

After  system  selection  and  optimization  was  accomplished  and  system 
calibration  completed,  a  data  run  was  conducted.  The  system  is  shown  in 
Figure  (14) .  The  data  run  was  commenced  with  the  variable  speed  trans- 
mission set  at  the  lower  value  of  rpm  corresponding  to  the  lower  value  of 
the  range  of  desired  frequency  parameters.  The  inductance  motor  and 
recorder  were  turned  on  simultaneously.  The  recorder  was  allowed  to  run 
until  the  channel  recording  wave  height  indicated  a  steady  state  had  been 
reached.  The  recorder  was  then  secured  and  the  transmission  adjusted  to 
the  next  rpm  value  corresponding  to  the  next  desired  frequency  parameter. 
The  meter  was  then  sevured  and  the  wave  action  in  the  channel  allowed  to 
subside  before  the  next  data  point  was  recorded.  This  procedure  was 
repeated  until  the  desired  range  of  frequency  parameters  had  been  covered. 

At  the  end  of  each  data  run,  a  calibration  check  was  conducted  on 
both  the  wave  height  gauge  and  strain  gauge  bridge. 
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IV.   EXPERIMENTAL  RESULTS 

As  shewn  in  Equation  (3)  in  THEORETICAL  CONSIDERATIONS,  the  force 
coefficient  can  be  represented  as  a  function  of  relative  water  depth, 
relative  depth  of  submergence,  relative  amplitude,  and  a  dimensionless 
frequency  parameter.  The  dynamic  forces  on  the  object  may  be  expressed 
as  the  sum  of  two  components,  one  component  in  phase  with  acceleration 
and  one  component  in  phase  with  the  velocity  of  the  cylinder.  The 
coefficient  of  the  acceleration  component  represents  added  mass  and  the 
coefficient  of  the  velocity  component  represents  damping.  It  follows 
that  these  two  coefficients  are  determinable  by  measuring  the  force  on 
the  cylinder  and  the  phase  angle  between  force  and  acceleration. 

Additionally,  if  a  control  volume  were  drawn  around  the  system,  it 
could  be  shown  through  a  consideration  of  conservation  of  energy  that 
the  wave  energy  flux  out  must  equal  the  energy  associated  with  damping 
if  viscous  effects  are  disregarded.  Therefore,  wave  height  is  directly 
related  to  the  degree  of  damping. 

Accordingly,  a  program  was  designed  to  measure  forces,  phase  angles 
between  force  and  acceleration,  and  wave  heights  within  a  range  of  inde- 
pendent dimensionless  variable  limited  by  time  and  equipment  available. 

Twelve  data  sets  were  taken  at  six  variations  of  relative  water  depth 
(d)  and  relative  depth  of  submergence  (h)  .  The  six  variations  were: 
d  =  2.0,  h  =  4.0;  d  =  2.0,  h  =  5.0;  d  =  3.0,  h  =  5.0;  d  =  3.0,  h  =  6.0; 
d  =  4.0,  h  =  6.0;  d  =  4.0,  h  =  7.0.  All  data  sets  were  taken  at  a 
relative  displacement  (x)  of  0.1.  The  range  of  dimensionless  frequency 
parameters  (v)  for  each  run  was  from  0.04  to  1.7.  The  low  frequency  half 
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of  each  run  was  conducted  with  a  six  inch  OD  plexiglass  cylinder  and  the 
high  frequency  half  with  a  four  inch  OD  aluminum  cylinder.  An  overlap 
between  high  frequency  and  low  frequency  parameters  of  approximately  0.36 
was  provided  near  mid-range. 

A.  EXAMPLE  STRIP  CHART  READ-OUT 

An  example  of  a  strip  chart  read-out  is  presented  in  Figure  (15) . 
For  this  particular  run  a=  3,  d=  6.0,  h  =  12.0.  The  values  read  frcm 
the  chart  for  this  data  point  were  wave  height  (WH)  =  0.088  inches, 
force  =  0.3525  pounds,  period  =  143  millimeters,  phase  angle  (PA)  =  8 
millimeters,  and  chart  speed  =  100  millimeters  per  second. 

The  wave  profiles  are  nearly  sinusoidal.  The  natural  frequency  of  the 
entire  mechanism  is  seen  superimposed  upon  the  force  profile.  It  is 
apparent  that  the  natural  frequency  is  considerably  greater  than  the 
forcing  frequency.  Difficulty  was  encountered  in  reading  phase  angle 
from  the  chart.  A  chart  speed  of  100  millimeters  per  second  was  required 
to  resolve  phase  angle.  Because  of  the  superposition  of  natural  fre- 
quency upon  the  force  profile,  it  was  necessary  to  "fair"  the  force 
profile  to  properly  read  phase  angle.  The  result  of  reading  phase  angle 
in  this  manner  led  to  inaccuracies,  especially  at  lower  frequencies  where 
forces  were  small  and  resulting  force  profiles  relatively  flat. 

B.  DATA  REDUCTION 

A  computer  program  was  written  to  reduce  experimental  data.  The 
program  is  presented  in  Appendix  B  with  tabulated  results  for  each  data 
point. 
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C.   PRESENTATION  OF  RESULTS 

Experimental  results  are  presented  in  graphical  form  in  Figures  (16-23) . 
For  each  combination  of  the  dimensionless  parameters  d  and  h,  a  graph  of 
force  coefficient,  phase  angle,  and  wave  height  ratio  was  drawn  to  exhibit 
dependence  upon  the  frequency  parameter.  In  a  few  cases  theoretical  curves, 
provided  by  Garrison  in  an  unpublished  paper,  were  drawn  in  for  comparison. 
Smooth  curves  were  "faired"  in  cases  where  Garrison's  theoretical  data  was 

not  complete. 

1.  Force  Coefficient 

Despite  the  scatter  in  the  values  of  force  coefficient,  results 
compared  favorably  with  Garrison's  theoretical  curves,  except  at  the  lower 
end  of  the  frequency  spectrum.  See  Figures  (18) ,  (22) ,  and  (23) . 

A  consideration  of  the  boundary  value  problem  associated  with  the 
problem  resulted  in  the  conclusion  that  the  free  surface  boundary  would 
act  as  a  solid  boundary  at  very  low  frequencies  and  as  a  constant  pressure 
boundary  at  very  high  frequencies.  It  follows  that  force  coefficient  as 
a  function  of  frequency  should  approach  a  constant  value  asymtotically  at 
both  the  high  frequency  and  low  frequency  end  of  the  data  range  and  vary 
in  a  continuous  manner  between.  Also,  since  the  constant  pressure  boundary 
at  the  high  frequency  end  tends  to  eliminate  pressure  differences  in  the 
fluid,  it  follows  that  the  constant  asymtotic  value  at  the  higher  frequency 
end  should  be  lower.  Further,  it  follows  that  as  the  relative  proximity 
of  the  free  surface  increases,  the  effect  should  be  more  pronounced.  The 
results  as  presented  in  Figures  (16,  18,  20-23)  indicate  that  these  effects 
were  observed. 
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Figure  (16)  indicates  that  for  d  =  2.0  and  h  =  4.0  the  lower 
frequency  value  of  force  coefficient  was  approximately  1.5.  The  range 
of  frequency  parameters  did  not  extend  to  high  enough  values  to  deter- 
mine the  high  frequency  asymptote.  Theoretical  data  was  not  available 
for  this  figure,  therefore  a  mean  line  was  drawn  through  the  data.  Force 
coefficient  appeared  to  reach  a  minimum  value  of  0.8  at  a  frequency 
parameter  of  1.15. 

At  the  greatest  depth  of  submergence/  (d  =  4.0,  h  =  6.0) ,  Figure 
(22)  indicates  that  the  lower  frequency  value  of  force  coefficient 
decreased  to  1.38.  Again,  the  high  frequency  asymptotic  behavior  was 
undeterminable.  Garrison's  theoretical  data  at  these  parameters  indicate 
that  the  force  coefficient  is  still  lower  than  experimentally  indicated 
throughout  almost  the  entire  range  of  frequency  parameters  by  as  much  as 
10%.  The  minimum  value  of  force  coefficient  at  this  larger  submergence 
depth  was  approximately  1.0  at  a  frequency  parameter  of  1.10. 

The  effect  of  the  proximity  of  the  impermeable  bottom  surface 
can  be  observed  by  comparing  Figures  (16)  and  (18) .  During  these  two 
data  runs,  relative  submergence  depth  (d)  was  held  constant  at  d  =  2.0 
while  relative  water  depth  (h)  was  varied.  Similar  comparisons  were 
made  at  d  =  3.0  and  d  =  4.0.  In  each  case,  the  proximity  of  the  impermeable 
bottom  surface  appeared  to  increase  the  force  coefficient  over  the  entire 
range  tested. 

It  may  be  noted  that  the  scatter  in  values  of  force  coefficient 
was  greatest  at  the  low  frequency  end  of  the  range  of  frequency  coeffici- 
ents. This  was  attributable  to  the  nature  of  the  error  in  force  measure- 
ment. It  was  observed  that  although  the  error  in  force  measurement  could 
be  assumed  relatively  constant,  the  magnitude  of  force  measured  at  values 
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of  v<  0.1  was  very  small,  often  smaller  by  an  order  of  magnitude  than 
the  magnitude  of  the  force  measured  at  values  of  p>  0.1.  It  follows 
that  with  constant  absolute  error,  the  relative  error  can  be  expected  to 
be  an  order  of  magnitude  larger.  In  every  case  scatter  was  reduced  at 
higher  frequencies  where  forces  were  larger.  Also  in  the  area  where 
overlap  occurred  between  high  and  low  frequency  data  sets,  the  results 
appear  consistant  and  repeatable.  It  was  observed  that  for  values  of 
0.1,  almost  all  of  the  data  was  consistant  within  10%. 
2.  Phase  Angle  and  Wave  Height 

Values  of  phase  angle  and  wave  height  compared  favorably  with 
Garrison's  theoretical  curves.  See  Figures  (18),  (19),  (22),  and  (23). 
As  stated  earlier,  the  added  mass  coefficient  is  determinable  if  the  force 
coefficient  and  the  phase  angle  are  known.  Further,  a  theoretical  con- 
sideration resulted  in  the  conclusion  that  added  mass  coefficient  may  be 
computed  by  multiplying  force  coefficient  by  the  cosine  of  the  phase  angle. 
The  balance  of  the  force  coefficient  is  attributable  to  damping  and  con- 
sequently, the  damping  coefficient  may  be  obtained  by  multiplying  the 
force  coefficient  by  the  sine  of  the  phase  angle.  Also,  as  stated  earlier, 
from  a  consideration  of  conservation  of  energy,  wave  height  may  be  used 
as  a  measure  of  the  degree  of  damping. 

The  value  of  phase  angle  and  wave  height  at  the  extremes  of  the 
frequency  spectrum  was  predicted  by  intuitively  assuming  that  at  very  low 
frequencies  the  force  and  motion  would  be  in  phase,  and  at  very  high 
frequencies  damping  would  be  non-existant .  It  was  predicted  that  phase 
angle  and  wave  height  would  be  non-existant  where  force  coefficient 
approached  a  constant  value  asymptotically.  It  followed  that  force 
coefficient  and  added  mass  coefficient  would  be  synonymous  near  the 
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asymptotes.  The  effect  of  a  free  surface  was  predicted  to  be  more  pre- 
dominant at  lower  values  of  relative  depth  of  submergence.  The  results 
as  presented  in  Figures  (16-23)  indicate  that  these  effects  were 
observed. 

Figure  (16)  indicates  that  for  d  =  2.0  and  h  =  4.0  the  phase  angle 
was  zero  for  frequency  parameters  greater  than  1.5.  The  lower  frequency 
data  was  not  complete  for  phase.  However,  Figure  (17)  indicates  that  wave 
height  was  zero  for  frequency  parameters  less  than  0.4.  Figures  (16)  and 
(17)  show  phase  and  wave  height  reached  maximum  values  of  0.48  radians 
and  0.74  inches  respectively. 

At  the  greatest  depth  of  submergence/  Figure  (22)  indicates  that 
for  d  =  4.0  and  h  =  6.0  the  phase  angle  was  zero  for  frequency  parameters 
greater  than  0.9.  Again,  the  phase  angle  did  not  vanish  at  lower  fre- 
quencies but  wave  height  vanished  for  frequency  parameters  less  than  0.4. 
Maximum  values  were  about  0.15  radians  and  0.24  inches  for  phase  angle 
and  waveheight  respectively.  As  predicted,  values  were  considerably  less 
than  those  observed  at  a  smaller  relative  submerged  depth. 

The  effect  of  the  proximity  of  the  impermeable  bottom  on  the  phase 
angle  can  be  observed  by  comparing  Figures  (16)  and  (18) .  The  maximum 
value  appears  to  be  reduced  and  shifted  slightly  to  the  left.  A  similar 
reduction  and  shift  to  the  left  were  observed  in  the  maximum  value  of 
waveheight  between  Figures  (17)  and  (19) . 

Again,  despite  the  scatter  in  the  values  of  phase  angle,  the 
results  compare  favorably  with  Garrison's  theoretical  curves.  In  cases 
where  theoretical  curves  were  not  available  the  trend  in  data  was  "faired 
in". 
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All  the  ronarks  above  pertaining  to  measuring  small  magnitude 
forces  with  constant  absolute  error  apply  also  to  the  measurement  of 
phase.  Accurate  measurement  of  phase  angle  was  especially  difficult  in 
the  lower  frequency  ranges  where  force  profiles  were  flat.  The  exact 
point  where  the  force  changed  direction  was  difficult  to  determine 
accurately.  Consequently ,•  more  scatter  was  observed  in  the  phase  data. 
It  was  for  this  reason  that  no  attempt  was  made  to  plot  the  added  mass 
dependency  upon  frequency  parameter.  Rather,  the  resulting  sets  of  force 
coefficient  data  and  phase  data  were  considered  separately.  Added  mass 
is  presented  in  tabular  form  as  computer  program  output  in  Appendix  B. 
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V.  CONCXUSIQNS 

An  experimental  study  of  the  added  mass  of  a  right  circular  cylinder 
under  a  free  surface  has  been  conducted.  From  the  results  of  this  study, 
the  following  conclusions  were  drawn. 

1.  Results  agreed  within  10%  with  theoretical  results  proposed  by 
Garrison  in  an  unpublished  paper. 

2.  The  relative  proximity  of  a  free  surface  induced  a  pronounced 
variation  in  added  mass  coefficient  from  the  classical  value  of  1.0. 

3.  When  the  cylinder  was  oscillated  at  very  low  frequencies  and  at 
at  very  high  frequencies,  the  value  of  the  added  mass  coefficient  was 
constant.  This  constant  value  was  larger  at  a  low  frequency  than  at 
high  frequency  due  to  the  fact  that  at  low  frequency  the  free  surface 
acts  as  a  rigid  boundary. 

4.  The  proximity  of  an  impermeable  bottom  surface  had  the  effect 
on  the  value  of  the  added  mass  coefficient  of  raising  its  magnitude  at 
all  frequencies. 

5.  A  phase  angle  between  induced  force  on  the  cylinder  and  the 
acceleration  experienced  by  the  cylinder  was  observed  to  be  a  function 
of  frequency.  This  phase  angle  vanished  at  high  and  at  lav  frequencies. 

6.  The  proximity  of  either  a  free  surface  or  an  impermeable  bottom 
surface  had  an  effect  of  raising  the  magnitude  of  the  maximum  value  of 
phase  angle  and  shifting  the  occurrence  of  this  maximum  to  a  lower 
frequency  value. 
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7.  When  the  cylinder  was  oscillated,  induced  waves  were  observed. 
Wave  height  was  a  function  of  frequency  of  oscillation.  Induced  waves 
vanished  at  low  and  at  high  frequencies. 

8.  The  proximity  of  either  a  free  surface  or  an  impermeable  bottom 
surface  had  the  effect  of  raising  the  magnitude  of  the  maxiinum  value  of 
wave  height  and  shifting  the  occurrence  of  this  maxiinum  to  a  lower 
frequency  value. 
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APPENDIX  B      COMPUTER  PRORAM 

C  THIS  PROGAM  WAS  USED  TO  REDUCE  EXPERIMENTAL  DATA, 

C  DEFINITION  OF  PARAMETERS  USED  ARE  AS  FOLLOWS: 

C  T=PERIOD  OF  FORCING  FUNCTION( SECONDS ) 

C  SIGMA=ANGULAR  FREQUENCY  OF  FORCING  FUNCTION 

C  (RADIANS/SECOND) 

C  G-GRAVITATIONM  CONSTANT ! FEET/SECOND**2 ) 

C  HB=WATER  DEPTH ( I NCHES ) 

C  AB=CYLINDER  RADIUS! INCHESI 

C  RHO=WATER  DENS ITY! SLUG S/FEET**3 ) 

C  LB=CYLINDER  LENGTH (  INCHES ) 

C  DB=DEPTH  OF  SUBMERGENCE! INCHES ) 

C  WH=TROUGH  TO  CREST  WAVE  HE IGHT ( INCHES ) 

C  F=MAXIMUM  FORCE! POUNDS ) 

C  XB=HALF  AMPLITUDE  OF  MOT  ION! I NCHES ) 

C  TMM=PERIOD  READ  FROM  CHART ! MILLIMETERS ) 

C  PMM=PHASE  SHIFT  READ  FROM  CHART! MI LL IMETERS) 

C  CSPEED=CHART  S PE ED ( MI LL I  METERS/SECOND) 

C  PHASE=PHASE  SH I  FT! RADI ANS ) 

C  NCARD=NUMBER  OF  DATA  CARDS 

C  ANU=DIMENSIONLESS  FREQUENCY  PARAMETER 

C  WHR=DIMENSIONLESS  WAVE  HEIGHT  RATIO 

C  FC=DIMENSIO:MLESS  FORCE  COEFFICIENT 

C  AM=ADDED  MASS  COEFFICIENT 

C  DN=DAMPING  COEFFICIENT 
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APPENDIX  B  (CONTINUED) 

10  FORMAT(4F10.4,  110) 

15  FORMAT ( •  1« ,20X, • AB  = • , F 8 .4, 5X , « DB= • , F8.4, 5X , »HB=f tF8.4, 

15X,,XB=,,F8.4,5X,1 NCARD=« , 13) 
20  FORMAT ( ' 0' , 20X , • X= • , F8 .4 , 5X , « D=« , F8.4, 5X , * H= • , F8 .4) 
25  FORMAT(5F10.4) 
65  FORM  AT (  ■ 0'  , 7X,  ' S IGMA'  , 10X , • F • , 7X , • C SPEED1  » 8X t • ANU»  f  10X 

If  •FC,,6X,  'PHASE'  ,10X,  •  AMSSXf'DN'  ,9X,  »WHR»  ) 
80  FORMAT! '0« ,9F12.4) 

ASSIGN  SPECIFIC  VALUES  TO  CONSTANTS  TO  BE  LATER  USED 

G=32.174 
RH0=62.4/G 

THIS  PROGRAM  WAS  DESIGNED  TO  HANDLE  TWELVE  DATA  SETS 

DO  300  K=l,12 

READ    IN    PARAMETERS    WHICH    CHARACTERIZE    DATA    SET 

READ(5,10)     AB,DB,HBtXB,NCARD 

PRINT    PARAMETERS    FOR    LATER    REFERENCE 

WRITE(6,15)    AB,DB,HBf XB,NCARD 

COMPUTE    VALUES    OF    DIMENSIONLESS    PARAMETERS 

D=DB/AB 
H=HB/Ac 
X=XB/AB 

PRINT    DIMENSIONLESS    PARAMETERS 

WRITE(6,20)     XfO,H 

PRINT    COLUMN    HEADINGS    FOR    INPUT    /    OUTPUT    PARAMETERS 

WRITE{6,65) 

COMMENCE  DATA  REDUCTION 

DO  200  I=1,NCARD 

READ  INPUT  PARAMETERS 

READ(5,25)  F , WH , TMM, PMM, CSPEED 

COMPUTE  PERIOD  OF  FORCING  FUNCTION 

T=TMM/CSPEED 

COMPUTE  PHASE  ANGLE 

PHASE=2.*3.14159*PMM/TMM 

COMPUTE  ANGULAR  FREQUENCY  OF  FORCING  FUNCTION 

SIGMA=2.*3.14159/T 

COMPUTE  DIMENSIONLESS  FREQUENCY  PARAMETER 

ANU=SIGMA**2*AB/(G*12) 
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APPENDIX  B  (CONTINUED) 

C      COMPUTE  WAVE  HEIGHT  RATIO 

WHR=(WH/2. )/XB 
C      COMPUTE  VOLUME  OF  CYLINDER 

VOL=3.14159*AB**2*l. 25/144, 

C      COMPUTE  FORCE  COEFFICIENT 

FC=F/(RHO*VOL*SIGMA**2*XB/12. ) 
C      COMPUTE  ADDED  MASS  COEFFICIENT 

AM=FC*COS( PHASE) 
C      COMPUTE  DAMPING  COEFFICIENT 

DN=FC*SIN(PHASE) 

C  COMPUTATION  COMPLETE.  PRINT  VALUE  OF  PARAMETERS  AT 
C  THIS  POINT  AND  RETURN  TO  BEGINNING  OF  LOOP  TO  LOOK 
C      AT  THE  NEXT  DATA  POINT  IN  THIS  SET. 

WRITE (6, 80)  SIGMA, F , C SPEED , ANU , FC , PHASE , AM, DN, WHR 
200  CONTINUE 

C      GO  TO  NEXT  DATA  SET. 

300  CONTINUE 
STOP 
END 
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Experimental  data  was  reduced  using  a  computer  program  and  the  results  presented 
graphically  as  a  function  of  dimensionless  parameters. 
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